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Abstract 

We study B — > <j)K and B (pXg decays in the heavy quark limit using per- 
turbative QCD. The next leading order corrections introduce substantial mod- 
ifications to the naive factorization results (more than 50%). The branching 
ratio Br{B 4>K) is predicted to be in the range (-F]^^^(m^)/0. 33)^(3. 5 ~ 
4.2) X 10^^ which is within the one a allowed region from the central value of 
6.2 X 10^^ measured by CLEO, but outside the one a allowed region from the 
central value of 17.2 x 10~^ measured by BELLE for reasonable F^'^^ . For 
the semi-inclusive decay B — > (pXs we also include initial bound state effect 
in the heavy quark limit which decreases the branching ratio by about 10%. 
Br{B (pXs) is predicted to be in the range (5.1 ~ 6.3) x 10^^. 
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Recently CLEO and BELLE have measured the penguin induced AS* = 1 hadronic B 
decays with Br{B- cjyK) = (6.411:51^;^ x IQ-^), Br{B^ (f)K^) = (S.Qttglag x IQ-^) 
from CLEO and5r(5- ^ (l)K') = {1.72toiltoilxlO~^) fromBELLE §. Although the 
central values do not agree with each other, they are consistent at 2a level. The branching 
ratios will soon be determined with better precisions at CLEO, BABAR and BELLE. These 
decay modes are particularly interesting in the Standard Model that they are purely due to 
penguins to the leading order and therefore are sensitive to new physics at loop level p . 
The neutral decay mode also provides a model independent measurement for one of the KM 
unitarity triangle parameter sin 2(3. The related semi-inclusive decay mode B —>■ (pXg is also 
purely due to penguin and is sensitive to new physics at loop level. The branching ratio 
for this decay although not measured at present, it will be measured in the near future at 
B factories. The above exclusive and semi-inclusive decays have been studied theoretically 
before with large errors 0,^J^ that both the CLEO and BELLE measurements can be 
accommodated. 

Previous calculations for the branching ratios Br{B (t>K) and Br{B (pXg) are 
based on naive factorization calculations. In these calculations, non-factorization effects can 
not be calculated and are usually parameterized by an effective color number and treated 
as a free parameter. There are also uncertainties related to gluon virtuality in the penguin 
diagrams and dependence of the renormalization scale. To have a better understanding of 
these decays, it is necessary to carried out calculations in such a way that the problems 
mentioned and others potential problems can be dealt with. It has recently been shown that 
it is indeed possible in the heavy quark limit to handle most of the problems mentioned in 
related B to two light mesons from QCD calculations . Several decays have been studied 
with interesting results PJIOj]. In this paper we will follow the method developed in Ref. |Q] 



to carry out calculations for the branching ratios for B —>■ (j)K and B —>■ (pXg. 
The effective Hamiltonian for charmless B decays with AS = 1 is given by 

^^e// = -4 j KfeC(ciOi + C2O2 + E c^o,) + v,bv:, J2 c^o. \ . 



(1) 



Here Oi are quark and gluon operators and are given by 

Oi = {SaUi3)v-A{upba)v-A, O2 = (SaU^) V-a(M/3&/3) V-A, 

03(5) = {Saba)v-A^{q'/3q'f3)v~{+)A, O^q) = (s„6/3) V-A 

q' q' 

3 3 
'^m = 2^Saba)v-A^ eg>{qpq'p)v+{-)A, Os{io) = -{sabf3)v~A^eqi{qpqJv+(-)A, 
q' q' 

On = ^,m,s^o^^G%^{l + ^,)bp, (2) 

where {V ± A){y ± A) = 7'^(1 ± 75)7^(1 ± 75), q' = u,d,s,c,b, tqi is the electric charge 
number of g' quark, Aa is the color SU(3) Gell-Mann matrix, a and /3 are color indices, and 
Gy,u is the gluon field strength. 

The coefficients q are the Wilson Coefficients which have been calculated in different 
schemes P,pT^. In this paper we will use consistently the NDR scheme. The values of q at 
/i ^ GeV with the next-to-leading order (NLO) QCD corrections are given by [jll 



ci = -0.185, C2 = 1.082, C3 = 0.014, C4 = -0.035, C5 = 0.009, ce = -0.041, 

cy = -0.002/aem, cs = 0.054/aem, cg = -1.292/aem, Clo = -0.263/aem, Cii = -0.143. 

Here a^m = 1/137 is the electromagnetic fine structure constant. 
The exclusive B ^ (f)K decay 

In the heavy quark limit, the decay amplitude due to a particular operator can be 
represented in the form ^ 

< <f)K\0\B >=< <PK\0\B >fact [1 + E ^na: + O(AqcdM)], (3) 

where < (j)K\0\B >fact indicates the naive factorization result. The parameter Aqcd ~ 
0.3GeV is the strong interaction scale. The second and third terms in the square bracket 
indicate higher order and h^Qco/^b corrections to the matrix elements. Including the 



next-leading order corrections and use information from Ref. []10[, we have the decay ampli- 
tude for B — > (j)K in the heavy quark limit 
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v2 

VcbV:,[al + al + al- ^{a'j + + a^o) + a^oa)]- 



(4) 



We will use the notation < 0|s7^6|i? >= m^f^ef^ and < K|s7^6|i? >= Ff^^ {q'^){p'^ + 

The coefficients a"''^ are given by 

a,=a, = c, + - + -—c,F^, 

al = Ci + ^ + [c3{F^ + G^{ss) + G^{sb)) + ciG^{sp} 

iV 47r iV 

b 

+ (c4 + ce) ^ G^{sf) + ciiG^,!! 
/=« 

«5=«5 = C5 + ^ + ^^Ca(-F,-12), 

a^ = 4 = c. + | + ^^cs(-F,-12), 

« c Clo 
a9=a9 = C9 + — + -— cioF^, 

u _ c _ Cg asGp „ 

Oio - ^10 - + ]7 + 4^ aT^^ 



"lOa 



as Gf 



4tt N 



3 ^ 3 
(c8 + cio)- ^ efG^{sf) + cg-{esG^{ss) + ebG^{sb)) 

f=u 



(5) 



where p takes the values u and c, = 3 is the number of color, Gp = (iV^ — 1)/2N, and 
Sf = m'j/ml. The other items are given by 



G^{i 
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/o 1 — X 
F. = -121n^-18 + /i + /, 
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Here are meson wave functions. In this paper we will take the following forms for them 
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(l>K,4,{x) = 6x{l - x), (7) 

where A^^ is a normalization factor satisfying dx(j)B{x) = 1. Fitting various B decay data, 
ujb is determined to be 0.4 GeV. 

The above results are from genuine leading QCD calculation in the heavy quark limit. 
The number of color should not be treated as an effective number, but has to be 3 from QCD. 
The results are renormalization scale independent. The problem associated with the gluon 
virtuality k"^ = (1 — a;)m^ in the naive factorization calculation is also meaningfully treated 
by convoluting the x-dependence with the meson wave functions in the functions G{s,x). 
Also leading non-factorizable is included (by the term proportional to f'l')- There are still 
uncertainties in the calculation, such as the form of the wave functions and the unknown 
B ^ K transition form factor Ff^^lq"^). However using wave functions obtained by fitting 
other data, the errors can be reduced. In any case calculations based on the method used 
here is on more solid ground compared with previous calculations. 

The decay rate can be easily obtained and is given by 



r(5 - <PK) = ^\C\^fiFr^{mirmlxZ (8) 



where Ajj = (1 — mf/m\ — m?-/m?^Y — Amfm'j/'m'^. 

In our numerical calculations we will use the following values for the relevant parameters 



T^: mb = 4.8 GeV, = 1.4 GeV, = 0.2196, Kb = 0.0395, Kfe/Kf, = 0.085, = 0.233 
GeV, fx = 0.158 GeV, and /b = (180 ± 20) MeV. We keep the phase 7 to be a free 
parameter. The results on the branching ratios are not sensitive to light quark masses. We 
obtain the branching ratios for B — > (pK to be 

BriB- <PK-) = (^^^^^^)' (3.7 ~ 4.2) x 10-^ 
Br{B' <f>K') = I ^^gg 1 (3.5 ~ 4.1) x 10"^ (9) 



We have checked sensitivities on some of the parameters. The branching ratios are 
insensitive to the phase angle 7 because terms proportional to e~*^ are suppressed by 
|\4fe\4*g/V^f,V^* I which is about 1/50. The error associated with Vcb is about 5%. The ranges 
given for the branching ratios have taken the dependence on 7 and the uncertainty in Vcb 
into account. We find that the NLO corrections to the matrix elements (terms proportional 
to as in tti) to be significant. Without such NLO corrections, the branching ratios are in the 
range of {F^~*^ {m'^)/ 0.33)"^ {2.3 ~ 2.5) x 10"^. The non-factorizable contributions (terms 
proportional to /^^) tend to reduce the branching ratios at a few percent level. The from 
factor Ff^^ is the least known parameter in the calculations. There are several calculations 



for this parameter. Lattice calculation gives 0.27 ± 0.11 [|T3|, BSW model gives 0.38 [jT4 



while light-cone calculation gives 0.35 ±0.05 [T^. Using the average central value from these 
calculations, Ff^^(O) = 0.33, one finds that the predicted branching ratios are closer to 
the averaged central value of the measurements from CLEO than that from BELLE. To 
reach the CLEO central values, Ff^^ needs to be around 0.42 which is on the high value 
side from theoretical calculations, while to reach BELLE central value a unreasonably large 
value 0.72 for Ff^^ is needed. Precise measurements of these modes may provide a good 
measurement of the form factor Ff~*^. If a better understanding of the form factor Ff~*^ 
can be obtained from other experimental measurements and from theoretical calculations 
in the future, precise measurement of i? ^ (pK may provide us with important information 
about new physics beyond the SM. 
The semi-inclusive B — > (pXs decay 

We will follow the procedures for semi-inclusive B decays described in Ref. |jl6|] to study 
B — > (f)Xs. The final state Xg can be viewed as containing a perturbatively produced s quark 
and some non-perturbatively produced state X containing no strange number. Neglecting 
color octet contribution the decay width for each of the helicity state A of 0(A), at the 
leading order with light quark masses set to zero, can be written as 



• < B\h^{l - 75)7 ■ qiuil - l5)Kx)\B > (10) 

• < 0|s(0)7^(l - 75)s(O)|0(A) + X X 0(A) + X\s{x)Yil - 75)s(x)|0 >, 

X 

where the parameter C is related to the Wilson Coefficients q. In the vacuum saturation 
approximation, 

J2 < 0|s(0)7'^(l - 75)s(O)|0(A) + X>< 0(A) + X\s{x)r{l - l5)s{x)\0 > 



X 



0|s(0)7'^(l - 75)s(O)|0(A) X 0(A)|s(x)7^(l - ^,)s{x)\0 > . (11) 

In this approximation the color octet contributions are automatically neglected. We 
will work with this approximation to estimate the branching ratio for B (pXg. This 
approximation is consistent with the assumption made in the previous section if color octet 
is neglected. If one cuts the momentum to be above 2 GeV or so, the contributions are 
dominated by the effective two body decay b 0s. In this case C is similar to C but 
with /^^ set to be zero. In principle terms proportional to /^^ also contribute. However 
this contribution is small and can be neglected. This is because that in the semi-inclusive 
decay only in the final state is specified. When the constraint of having K in the final 
state is relaxed, the term corresponding to /^^ leads to a three body decay. Requiring the 
identified hadron in the final state to be hard limits the phase space |T6| and results in a 
small contribution from /^^ compared with other contributions. 

If the b quark mass is infinitively large, Br{B (pXs) is equal to Br{b — > 0s). However 
due to initial b quark bound state effect there are corrections . This correction is included 
in the factor < B\b{0)j^{l — 75)7 ■ ^7*^(1 — 75)6(x)|i? >. Following the discussions in Ref. 



T7| we obtain 1 / mf, correction factor 



where 



TiB -> 0X.) ^ ^^mi + 1^- f 4)> (12) 



/ij = < B\h^g,G^,a^^''h\B >, 

l4 = - < B\hDlh\B > . (13) 



Here the field h is related to b by b(x) = e"*"*'^'^ {1 + ^7 ■ DT/2mb + v ■ ■ DT/^ml — (7 • 
DtY /Sml}h{x) + 0(l/m^) + (terms for anti-quark). Dj, = — v^v ■ D with v being the 
four velocity of the B meson satisfying f ^ = 1 and = + igsG^{x). 



Analysis of spectroscopy of heavy hadrons and QCD sum rule calculations give |]T8 
^ 0.36 GeV^ and fil ^ (0.3 ~ 0.54) GeV^. We will use ^^ = 121 = 0.36 GeV^ for 

numerical calculations. The initial state 1/ml correction reduces the branching ratio by 

about 10%. The branching ratio for B (pXg is predicted to be 

Br{B (f)Xs) = (5.1 ~ 6.3) x 10"^ (14) 

This prediction, as in the case for exclusive decays, is insensitive the phase angle 7. The 
NLO corrections enhance the branching ratio significantly, similar to the exclusive decay 
cases. 

The expression for the semi- inclusive decay in eq.(0), on the face of it, has less param- 
eters (no dependence on Ff~^^) compared with the exclusive branching ratios discussed 
earlier. One might think that the prediction for Br{B ^X^) is more certain compared 
with the exclusive cases. However, one should be careful about this because in the calcu- 
lation we have only included color singlet and the ^Si ss bound state contribution. There 
may be other contributions such as color octet and other S and P wave states from ss. 
These contributions are in general smaller than the contributions already considered. One 
can not rule out significant enhancement at the present. However, we can view the color 
singlet result as leading contribution which gives a good order of magnitude estimate of the 
semi-inclusive decay B — » 4>Xs. 

In conclusion, we have studied B (pK and B (pXg decays in the heavy quark 
limit using perturbative QCD. We found that the next leading order corrections introduce 
substantial modifications to the leading native factorization results (more than 50%). The 
branching ratio Br{B (j)K) is predicted to be in the range (-F]^^^(m0)/0. 33)^(3. 5 ~ 
4.2) X 10^^ which is within the one a allowed region from the central value of 6.2 x 10~^ 
measured by CLEO, but outside the one a allowed region from the central value of 17.2 x 10~^ 
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measured by BELLE for reasonable Ff^^. For the semi-inclusive decay B — > (f)Xs we also 
included initial bound state effect in the heavy quark limit which decreases the branching 
ratio by about 10%. Br{B — > (pXg) is predicted to be in the range (5.1 ~ 6.3) x 10~^. 
Future experimental data will provide us with more information about these decays and 
about method based on QCD improved factorization approximation. 

This work was supported in part by NSC under grant number NSC 89-2112-M-002-016. 
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